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Damping of giant resonances in asymmetric nuclear matter 

00 ; 

K. Morawetz a,b , U. Fuhrmann a and R. Walke a,b 

p ■ a Fachbereich Physik, Universitat Rostock, 18051 Rostock, Germany 
D ■ b Laboratorio Nazionale Del Sud, Via S. Sofia, 44-95123 Catania, Italy 

m ; 

The giant collective modes in asymmetric nuclear matter are investigated within a 

dynamic relaxation time approximation. We derive a coupled dispersion relation and 

show that two sources of coupling appear: (i) a coupling of isoscalar and isovector modes 

due to different mean-fields acting and (ii) an explicit new coupling in asymmetric matter 

due to collisional interaction. We show that the latter one is responsible for a new mode 

Q\ ' arising besides isovector and isoscalar modes. 
O ! 
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1. Theoretical Preliminaries 

43 : 

| The treatment of collective modes in nuclear matter is documented in an enormous 

literature [ |J. We want to investigate the collective excitations in asymmetric nuclear 
matter [ 0, §. To this end we like to focus on the simplest microscopic theory which 
provides basic experimental features. This will be accomplished by a Fermi-gas model 
including dissipation. We consider the system consisting of a number of different species 
(neutrons, protons, etc..) interacting with the own specie and with other ones. It is 
important to consider the interaction between different sorts of particles if we want to 
include friction between different streams of isospin components, etc. Especially the 
isospin current may not be conserved by this way. Let us start with a set of coupled 
quantum kinetic equations for the reduced density operator p a for a specie a 
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d tPa (t) = i[ Pa , £ a + U{t) a \ - W dt' Pa{t ') y (1) 

b J Q T ab{t t ) 

where S = V 2 /2m denotes the kinetic energy operator and IA the mean field operator 
and the external field which is assumed to be a nonlinear function of the density. We 
approximate the collision integral by a non-markovian relaxation time. This turned out 
to be necessary to reproduce damping of zero sound [ [|, |J. It accounts for the fact 
that during a two-particle collision a collective mode can couple to the scattering process. 
Consequently, the dynamical relaxation time represents the physical content of a hidden 
three-particle process and is an equivalent to the memory effects. Further we assumed a 
relaxation with respect to a local equilibrium which we specify by a small deviation of the 
chemical potential of specie a determined by current conservation. Linearizing the kinetic 
equation (fj) we obtain the matrix equation for the density deviation n b with respect to 
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an external perturbation £7g Xt 
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The matrix C a b is given by 
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with a a b = ~^U a the linearization of the mean field with respect to the deviation of 



density from equilibrium value caused by the external perturbation U ex t- 
polarization function of a specie a is 



dp fa{P+l)-fa{P 
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The partial 
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The factor 2 in front of the integral accounts for the spin degeneracy according to n a = 

2/^/a(p). 

Equation (§) represents the complete polarization of the system, because n = UU ext - It 
represents a matrix equation which is solved easily. The collective modes are given by the 
zeros of the determinant of the matrix on the left hand side of @ because these are the 
poles of the polarization matrix. Since we took into account relaxation processes between 
all species we are able to cover current - current friction. For a two component system, 
e.g. neutrons with density n and protons with density p, we write explicitly 



[I 



C« nn )(i 



(,D n p <y n p^(Dp n -\- rip oip n j 



with the generalization of the Mermin polarization function 
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and the additional coupling due to asymmetry in the system 
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7| to a multicomponent 
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The D pn are given by interchanging sort indices. This term does not appear for symmetric 
matter. Therefore we call this term asymmetry coupling term further-on. 

It is illustrative to see known results for symmetric nuclear matter. This is performed 
for the case of equal relaxation times r„ = r n = r and equal deviation from the mean field 
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a pp and a2 
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relation [ 0] is then obtained 
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From (ft) the known Mermin result of dispersion 
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with the isovector mode oc\ — «2 and the isoscalar mode ot\ + a-i- 

We have presented a general dispersion relation for the multicomponent system in- 
cluding all known special cases. The dispersion relation (||) is similar to the one derived 
recently in [ |8| if we neglect the collisional coupling D np . Also a more general polarization 
function @ is presented here including collisions within a conserving approximation [ |§ . 
In the following we will apply this expression for the damping of giant dipole resonances 
for symmetric as well as for asymmetric nuclear matter. 



2. Model for nuclear matter situation 



We connect the wave vector for giant resonances according to the the Steinwedel- Jensen 
[ |10| model with the nuclear radius. If we assume that the density oscillation obeys a 
wave equation we get for the boundary condition that the radial velocity vanishes on the 
spherical surface with radius R = 1.13A 1 / 3 such that j[{kR) = with the spherical Bessel 
function of order I = 0, 1... associated with the monopole, dipole... resonances. 

Since monopole modes are compression modes we have no zero of first order for / = 0. 
For the dipole modes one has in first order k = 2.08/i? which describes the giant isovector 
dipole resonance (IVGDR) while the ISGDR is a spurious mode in first order. This would 
just mean an unphysical oscillation of center of mass motion. However, in second order 
k = 5.94/ R the isoscalar giant dipole resonance (ISGDR) has been observed recently [ I I 
and citations therein. This can be considered as a density oscillation inside a sphere. The 
resulting wave vectors have very low values compared with the Fermi wave vector. This 
allows us to expand the Mermin polarization function (^) with respect to small qv c /u 
ratios and v c the sound velocity. We obtain 
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m a u{u + i/T a y 

n a (n a ) = 2\~ 3 f 3/2 (z a )Ccorr 

where the thermal wave length is X 2 = 2nh 3 / \m a T) , the standard Fermi integral 
and z a = e Ma//T the fugacity. The correction constant ( corr = 1.22 is introduced to fit the 
numerical solution of the dispersion relation with the approximative expansion (|S|). With 
the help of this expansion the dispersion relation (H) takes the form 
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with the partial sound velocities c and c 



C ab = a ab- 



m„ 



r 2 



/3/2(zq) _ h/2{zb) 
T /l/ 2 (Za) /i/a(«%) 



Tab T~ aa 



(12) 



The dynamic relaxation time has been derived via Sommerfeld expansion [ ^ as 
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for a, b neutrons or protons respectively. The markovian relaxation time was given in 
terms of the cross section <r ab between specie a and b as r^ 1 = ^a^T 2 . The dispersion 
relation (|5]) takes then the form of a polynomial of tenth (sixth) order corresponding to 
the inclusion of memory (in) dependent relaxation times. While most of these solutions 
will just lead to parasitary solutions (Re a; < 0) we will get two coupled modes, i.e. the 
isoscalar and isovector mode. Furthermore a third mode appears at extreme asymmetries 
and/or strong collisional coupling. 

The damping rate in classical approximation is given by the solution of @ with (||) 
(long wavelength) and reads 7 = l/(2r) . We recognize that the FWHM is just twice the 
damping rate V = 27. This has been recently emphasized [ |12 |. It has to be stressed that 
the experimental data are accessible by this FWHM. 

We will use as an illustrative example the following mean field parameterization of 
Vautherin [ [T3], |TJ 
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with a = n,p the density of neutrons or protons respectively. The Coulomb interaction 
leads to an additional contribution for the proton mean-field 

U°{q) = -^n p (q). (15) 

The here used model parameters reproduce the Weizsacker formula 

E CL2 (x^Z 2 9 / \ 

J = - a i + A^ + A^ + ai5 (16) 
by the volume energy ai = 15.68 MeV, Coulomb energy 03 = 0.717 MeV and the sym- 
metry energy 04 = 28.1 MeV with the asymmetry parameter 5 = " n ~ np . 



First we plot the solution of the dispersion relation (|1T|) for symmetric nuclear matter. 
In figure [I] we plotted the real and imaginary (FWHM) part of complex energy for different 
approximations with relaxation time (collisions) with and without Coulomb. In [figure [I] 
(A)] we find that the inclusion of Coulomb effects reproduces the experimental shape of 
the centroid energies at higher mass numbers (dot-dashed line). Taking only collisions into 
account fails to describe higher mass numbers (solid line). Considering Coulomb together 
with collisions (dashed line), the centroid energies are reduced towardes the data. 

The experimental values of the damping rates are also presented versus mass number 
[figure [IJ (B)]. Without collisions we would only have a vanishing Landau damping for the 
infinite matter model [ 0] . We see that the inclusion of Coulomb (dashed line) improves 
the situation. 

The failure of the considered processes becomes more drastic if we consider the temper- 
ature dependence of the GDR. In figure § we see that the temperature increase is too flat 
compared to the experimental finding if we consider only collisional damping (thin lines). 
In Ref. [[Il| we present a method to include also scattering with the nuclear surface. This 
improves the temperature dependence remarkably (thick lines). 

Isoscalar giant dipole resonances (ISGDR) are observed recently in 208 Pb [ |ll|] and are 



given with a value 

£^ GDR = 22.4 ± 0.5 MeV, rg GDR = 3.0 ± 0.5 MeV. (17) 
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Figure 1. The experimental centroid 
energies (A) and damping rates (B) 
[filled symbols] of the giant dipole reso- 
nances vs mass number (data from Ref. 
[ 15 1) together with different approxi- 
mations at T = 0. 
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Figure 2. The effective damping 
consisting of only collisional damping 
compared with collisional and surface 
damping [ 120 Sn (solid lines) and 208 Pb 
(dashed lines)] together with the ex- 
perimental data (filled symbols: Sn 
from Ref. [ [T(| and Pb from Ref. 

E3)- 



They can be considered as a density oscillation inside a sphere. We obtain with the model 
parameters of IVGDR and a corresponding second order in j[{kR) = 

£isgdr = 21.0 MeV, r ISG DR = 3.4 MeV (18) 

which is quite acceptable and is in the middle of different theoretical treatments ([ [11]] 
and citations therein). 

3. New collective mode 

Now we employ the same potential as in the last section but assume different neutron 
and proton densities. In figure [| we plot the isoscalar and isovector modes versus tem- 
perature for 48 Ca with a small asymmetry 5 = 0.2 as well as for 60 Ca with a asymmetry 
(5 = 0.33. The kinetic energy is linked to a temperature within the Fermi liquid model 
via Sommerfeld expansion 



E 3 
A ~ 5 £/ 



1 + 5)5/3 + 5 )5/3 5 2 / T \ 2 ( 1 + 5 )l/3 + (1 _ 5) l/3 




(19) 



This connection between temperature and excitation energy is only valid for a continuous 
Fermi liquid model. For the small nuclei, the concept of temperature is questionable. 
Some improvement one can obtain by the definition of temperature via the logarithmic 
derivative of the density of states [ |19| 
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Figure 3. The energy and damping of the IVGDR and ISGDR of A8 Ca (left) and eo Ca 
(right) vs excitation energy. Besides the isovector modes (circles) and the isoscalar modes 
(squares) a soft third mode appears (triangles). 



which provides E ex « \{E/A) in comparison with (|TJ]) for small temperatures. We use 
this temperature to demonstrate possible collective bulk features in an exploratory sense. 
Of course, the surface energy and shell effects cannot be neglected for realistic calculations. 

With increasing temperature the isovector and isoscalar energies decrease and vanish 
at a certain temperature. At these temperatures the damping becomes twofold because 
the damping of the spurious mode with negative energy becomes different from the phys- 
ical mode. We can consider this behavior of damping as a phase transition of isospin 
demixture. 

At the same time a very soft mode appears due to collisional coupling which is only 
present in asymmetric matter [ |20[. This mode is more pronounced in the next figure ^| 



for m Ca with an asymmetry of 5 = 0.33. 

We see that the isovector mode does not disappear but turns over into a flat decrease 
with increasing temperature. This behavior is coupled with the pronounced soft mode. 
In comparison with the more symmetric 48 Ca we see a different behavior of the damping 
where the isoscalar mode vanishes. Also the isovector mode appears unique and not two- 
folded. Now a clear transition of damping behavior for the isovector mode is recognizable 
which can be considered as a transition from zero to first sound damping. 



4. Summary 

To summarize we have observed that due to correlational coupling there can exist a 
new mode which appears besides isovector and isoscalar modes in asymmetric nuclear 
matter due to collisions. We suggest that this mode may be possible to observe as a soft 
collective excitation in asymmetric systems. The transition from zero sound damping to 
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first sound damping behavior should become possible to observe for isovector modes since 
they do not vanish at this transition temperature like in symmetric matter. 
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